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ABSTRACT 

We investigate the star formation history of the universe using FIREWORKS, a multiwavelength 
survey of the CDFS. We study the evolution of the specific star formation rate (sSFR) with redshift in 
different mass bins from z = to z ~ 3. We find that the sSFR increases with redshift for all masses. 
The logarithmic increase of the sSFR with redshift is nearly independent of mass, but this cannot 
yet be verified at the lowest-mass bins at z > 0.8, due to incompleteness. We convert the sSFRs to 
a dimcnsionlcss growth rate to facilitate a comparison with a semi-analytic galaxy formation model 
that was implemented on the Millennium Simulation. The model predicts that the growth rates and 
sSFRs increase similarly with redshift for all masses, consistent with the observations. However, we 
find that for all masses, the inferred observed growth rates increase more rapidly with redshift than 
the model predictions. We discuss several possible causes for this discrepancy, ranging from field- 
to-field variance, conversions to SFR, and shape of the IMF. We find that none of these can solve 
the discrepancy completely. We conclude that the models need to be adapted to produce the steep 
increase in growth rate between redshift z = and z = 1 . 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift 



1. INTRODUCTION 

To understand galaxy formation and evolution, it is 
essential to have a clear idea of how galaxies assemble 
their mass. Multiwavelength galaxy surveys of the high- 
redshift universe provide estimates of masses and star 
formation rates out to z ~ 6 (Lilly et al. 1996, Madau 
et al. 1996, and the compilations by Hopkins (2004) and 
Hopkins & Beacom (2006)). Recent studies of the evolu- 
tion of the SFR per unit mass (specific SFR, sSFR) have 
shown that this quantity is an increasing function with 
redshift out to z ~ 2 and that, at a given redshift, the 
most massive galaxies typically have the lowest specific 
SFRs (Juneau et al. 2005, Bauer et al. 2005, Perez- 
Gonzalez et al. 2005; 2008, Caputi et al. 2006, Papovich 
et al. 2006, Reddy et al. 2006, Noeske et al. 2007, Mar- 
tin et al. 2007, Zheng et al. 2007, Damen et al. 2009). 
Many of these studies support the idea of downsizing, 
where the locus of active star formation shifts from mas- 
sive galaxies to less massive galaxies with time (Cowie et 
al. 1996). 

It is interesting to compare these observations to model 
predictions. Using semi-analytical modeling (SAM) tech- 
niques, the formation of galaxies can be simulated within 
the standard cold dark matter (ACDM) cosmogony. 
SAMs have been relatively successful in reproducing nu- 
merous systematic properties of the observed local pop- 
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ulation at z <0.4, such as number densities, luminos- 
ity functions, mass functions, and SFRs (e.g. Madau, 
Pozzetti & Dickinson 1998, Cole et al. 2000, Fontana et 
al. 2006, Fontanot et al. 2009). Since the recipes in the 
SAMs are arguably tuned to obtain such a good match, it 
is an interesting exercise to extend the comparison with 
observations to higher redshift, leaving the parameters 
of the SAMs to their locally tuned values. Several stud- 
ies make the comparison for the high-redshift universe 
(Bower et al. 2006; Kitzbichler & White 2007; March- 
esini & van Dokkum 2007; Dave 2008; Genzel et al. 2008; 
Guo & White in preparation, Elbaz et al. 2007, Daddi 
et al. 2007a, Santini et al. 2009). 

Bower et al. (2006) find their model to be in good agree- 
ment with the observed mass function out to z ~ 2. How- 
ever, their models overpredict the star formation density 
by ~ 20% at low redshift (z < 0.4), whereas at z ~ 1 
they underpredict the observations by a similar factor. 
Refinements to the code with respect to heating by ac- 
tive galactic nuclei (AGN; Bower et al. 2008) have not 
yet helped to reduce the relatively small discrepancies. 
A similar effect, where the models underpredict star for- 
mation observed at high redshift, is also seen in hydro- 
dynamical simulations (both for multiwavelength (Dave 
2008) and spectroscopic (Genzel et al. 2008) observa- 
tions. Kitzbichler & White (2007) find an abundance 
of massive galaxies that overpredicts the observed abun- 
dance by a factor of ^2 at z = 2. Studies by Elbaz et al. 
(2007), Daddi et al. (2007a), and Santini et al. (2009), 
reveal how SAMs underpredict the observed star forma- 
tion rate by a factor 2-5 out to redshift z = 2. 
In this paper we explore the evolution of the star forma- 
tion history of galaxies from redshift z ~ 3 to 0, using 
FIREWORKS, a .Kg-band selected multiwavelength cat- 
alog of the CDFS, that includes deep MIPS data to allow 
derivation of SFRs. Galaxy growth can be conveniently 
quantified using the dimensionless growth rate, defined 
as sSFR *tu(z), where tn(z) is the Hubble time (Guo & 
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Fig. 1. — sSFR against mass, for the low redshift sample from the SDSS, and at mean redshifts, z = 1,2,3. Open circles are the 
FIREWORKS detections, arrows show 2cr upper limits of the sSFR. 



White 2008). It provides a direct and quantitative con- 
straint on the models. We first compare our results to 
other observational studies and then, using the growth 
rate, make a comparison with a SAM built on top of the 
Millennium Simulation. Throughout this paper we as- 
sume a ACDM cosmology with £7 m = 0.3, f^A = 0.7, and 
Hq = 70 km s _1 Mpc -1 . All magnitudes are given in the 
AB photometric system. Stellar masses are determined 
assuming a Kroupa (2001) initial mass function (IMF). 

2. DATA 

2.1. Observations and sample selection 

We use the FIREWORKS catalog for the GOODS- 
CDFS, generated by Wuyts et al. (2008). This catalog is 
based on high quality photometry ranging from the NUV 
to MIR. Wuyts et al. (2008) included deep space based 
optical imaging with HST using ACS (Giavalisco, Stcidcl 
& Macchetto 1996). They complemented this with opti- 
cal imaging obtained as part of the COMBO-17 (Wolf et 
al. 2003) and ESO DPS (Arnouts et al. 2007) surveys, in 
the re-reduced form of the GaBoDs consortium (Erben 
et al. 2005; Hildebrandt er al. 2006). NIR imaging was 
obtained with VLT/ISAAC (Vandame 2002, Vandame 
et al. in preparation). For mid-IR wavelength coverage 
they used deep Spitzer imaging (IRAC and MIPS) from 
Dickinson et al. (in preparation). 

A K s -selected catalog was constructed following the pro- 
cedures of Labbe et al. (2003) and contains the follow- 
ing bands: U 38 BVRI (WFI), £435^06*775*850 (ACS), 
JHK S (ISAAC), 3.6-8.0 /^m (IRAC) and 24 /im (MIPS). 
It has a 5 sigma depth in K s of ^24.3 and a total area 
of 138 arcmin 2 . For details on observations, source de- 
tection and astrometry we refer to Wuyts et al. (2008). 
Using the CDFS X-ray catalog of Giacconi et al. (2002), 
we excluded all X-ray detected sources from the sample 
as they are likely AGN. We further restricted the selec- 
tion to sources with a signal-to- noise higher than 10 in 
the i^-band, which results in a total sample size of 5,274 
sources. 

2.2. Redshifts, masses and star formation rates 

Wuyts et al. (2008) derived photometric redshifts us- 
ing EAZY (Brammer et al. 2008) which are in good 
agreement with the available spectroscopic redshifts. 
The normalized median absolute deviation (NMAD) of 
(zphot - z sp ec)/(l + Zspcc) is 0.031 over the whole red- 
shift range and NMAD = 0.071 at z > 1.5 (Wuyts et 
al. 2008). To obtain stellar masses, ages, dust attenua- 
tion and star formation rates, stellar population models 



were fitted to the spectral energy distribution (SED) out 
to 8 ^m. The SED modeling follows standard proce- 
dures, and is described by Forster Schreiber et al. (in 
preparation; see also e.g., Forster Schreiber et al. 2004). 
Bruzual & Chariot (2003) models were used with solar 
mctallicity, a Salpeter IMF and a Calzetti reddening law. 
Masses and SFRs were converted to a Kroupa IMF by 
subtracting 0.2 dex. The fits were performed with the 
HYPERZ program (Bolzonella, Miralles & Pello 2000), 
fixing the redshift to the photometric redshift derived 
with EAZY, or the spectroscopic redshift when available. 
Three star formation histories were fit and the best fit- 
ting model was used. The three star formation histories 
were: a single stellar population without dust, a con- 
stant star formation history and an exponentially declin- 
ing star formation history with a characteristic timescale 
of 300 Myr, the latter two with varying amounts of dust 
(0 < A v < 4). A full description of the SED fitting pro- 
cedure and extensive tests on the outcome can be found 
in Forster Schreiber et al. (in preparation). Rest-frame 
luminosities were derived by interpolating between ob- 
served bands using the best-fit templates as a guide (see 
Rudnick et al. (2003) for a detailed description of this 
technique and Taylor et al. (2009) for the IDL imple- 
mentation of the algorithm, dubbed 'InterRcst' 7 ). 
In addition to the SFRs derived from SED fitting, SFRs 
were determined independently using a combination of 
rest- frame UV and IR emission. Between 1 < z < 3, 
MIPS 24 ^im traces the rest-frame 6-12 ^m luminos- 
ity, which correlates broadly with the total IR luminos- 
ity. Using a wide range of templates from Dale & Hclou 
(2002), bolomctric luminosities were determined, adopt- 
ing the mean of log(L/#), following Wuyts et al. (2008) 
and Labbe et al. (in preparation). In other words, we 
apply a linear conversion from observed 24 ^m flux to 
Lih, using a coefficient which depends also on redshift. 
The SFR was determined assuming 

*/M Q yr- 1 = 1.1 x 10- 10 x (L m + 3.3 L 2sm )/L Q , (1) 

based on the conversion by Bell et al. (2003) converted 
to a Kroupa IMF (see also Wuyts et al. (2008), Labbe 
et al. (in preparation)). Unless stated otherwise, we use 
the MIPS derived SFRs in our analysis. The SFR de- 
termined from the SED fitting is used for comparison in 
later sections. 

We include data from the Sloan Digital Sky Survey 
(SDSS) to extend our redshift coverage to z = 0. SDSS 

7 http://www.strw.leidenuniv.nl/~ent/InterRest 
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TABLE 1 

75% COMPLETENESS LIMITS 



log(Af„/Afo) 


z 


9.0 -9.5 


0.5 


9.5 -10.0 


0.9 


10.0 -10.5 


1.5 


10.5 -11.0 


2.9 


>11.0 


2.9 



masses were determined by Kauffmann et al. (2003) us- 
ing spectra. Brinchmann et al. (2004) derived SFRs 
from emission lines. For details on the derivation of the 
masses and SFRs in the SDSS we refer to their papers. 
Figure [T]shows the specific star formation rate with mass 
in four different redshift regimes [z = 0,1,2,3). The 
first panel shows SDSS results, the others are based on 
the FIREWORKS data. Although the spread in sSFR 
is high at a given mass (see also Franx et al. 2008), 
two trends are clear. Galaxies of higher mass typically 
have lower sSFRs and at fixed mass the sSFR rises with 
redshift. Both trends have been observed before and con- 
firm earlier results (Bauer et al. 2005, Zheng et al. 2007, 
Perez-Gonzalez et al. 2008, Damen et al. 2009). We will 
discuss the evolution of the sSFR further in Section [3J 

2.3. Mass completeness 

We determine the mass completeness limit using 
our i4T s -band selection limit. We scale the masses 
of sources down to flux selection limit of 10 a 
(M sca i ed — M x 10/[S/N] Ks ) and determine the 
mass limit to which we can detect 75% of the sources at 
10 a in a narrow redshift bin. We find that our sample 
is 75 % complete for masses higher than 1O 11 M0 out to 
z ~ 3. The completeness limits are listed in Table [TJ In 
the remainder of the paper, all references to complete- 
ness are based on the 75% completeness limits. 

3. THE EVOLUTION OF THE SPECIFIC STAR 
FORMATION RATE 

In a previous paper (Damen et al. 2009) we stud- 
ied how the SFR changes with mass and redshift out to 
z ~ 2. We used the SIMPLE survey, which is based on 
NUV to NIR observations of the E-CDFS. With an area 
of 0.5 x 0.5 deg, the E-CDFS is ~6.5 times larger than 
the much deeper FIREWORKS area. Here we extend 
this analysis to higher redshift and lower masses, using 
the FIREWORKS survey of the CDFS. We determine 
the average SFR in five different mass bins and inves- 
tigate the evolution of the specific SFR with redshift. 
For sources with no significant MIPS flux, we did not 
use upper limits to determine the average, but included 
the measured fluxes. Our results are shown in Figure [2j 
Different colors represent different mass bins and dots 
denote where we suffer from incompleteness. Error bars 
represent the bootstrapped 68% confidence levels on the 
average sSFRs. Local SDSS data points (triangles) were 
added to the figure and lie on the same trend as our re- 
sults. We also compare our results to previous similar 
studies. Dotted lines show the results of Zheng et al. 
(2007) who used the COMBO- 17 survey. Dashed lines 
represent the work of Damen et al. (2009), based on the 
SIMPLE survey. 



TABLE 2 

SLOPE OF SSFR-2 RELATION 



log(M»/M ) 


n 


9.5 -10 


4.40± 0.31 


10 -10.5 


3.36±0.87 


10.5-11 


3.63± 0.33 


>11 


4.78± 0.37 



Those papers showed that sSFRs of massive galaxies are 
typically smaller than those of less massive galaxies and 
that the increasing trend of sSFR with redshift has a 
similar slope for all masses. For the lowest to average 
mass bins (9.0 < log(A//Ajf ) < 11) we confirm results 
from Zheng et al. (2007) and Damen et al. (2009) and 
extend them out to higher redshifts. However, the sS- 
FRs in the highest mass bin (M* > 10 11 M Q ), seem to be 
lower than both the results from SIMPLE and Zheng et 
al. 2007. This can be a matter of low number statistics, 
since the average number of galaxies per redshift bin for 
these masses is ~ 8. 

To quantify the slopes of the trends in each mass bin, 
we fitted the sSFR with (1 + z) n over the redshift range 
where it is complete with respect to mass. We used a 
bootstrapping technique to determine errors on the fits. 
The resulting values for the steepness of the slope, n, are 
listed in Table [2j the fits themselves are shown in Figure 
[2] by dash-dotted lines, in gray. 

It is interesting to regard these values in the context of 
downsizing. We can see that at fixed redshift, the sSFR 
decreases with increasing mass. This indicates that the 
most massive galaxies have already formed the bulk of 
their stars and that active star formation has shifted to 
galaxies that are less massive. However, we do not see a 
strong mass dependence in the decline of the sSFR with 
time. This confirms results from Zheng et al. (2007) 
and Damen et al. (2009), whereas it seems to disagree 
with the results of Juneau et al (2005), who find that, 
out to z — 1.5, the sSFRs of more massive galaxies de- 
crease faster in time than those of the lowest-mass galax- 
ies. Taking into account the error bars in their Figure 3, 
this disagreement could be subtle. 

Admittedly, the redshift range over which we determined 
n is small for the low-mass bins (see Figure [2]). In ad- 
dition, the highest mass bin (log(.M*/ Mq) > 11 suffers 
from low number statistics. We verified whether our re- 
sults would change with a larger number of sources by 
including the SIMPLE data, and find that the trend in 
Figure[2]at the high-mass end remains the same. Includ- 
ing the SIMPLE data does not alter the results in a way 
that evolution in the slope can be confirmed. To be more 
conclusive, these fits need to be determined over similar 
redshift ranges for all mass bins. Deeper observations, 
that solve the incompleteness issues, are needed to accu- 
rately make a statement about the possible evolution of 
n with mass. 

4. COMPARISON WITH MODEL PREDICTIONS 

Next, we compare our results with model predictions. 
We use the work of Guo & White (2008), who have 
used data from the Millennium Simulation to study the 
growth of galaxies through mergers and star formation in 
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Fig. 2. — sSFR versus redshift in different mass bins. Filled, connected circles are the FIREWORKS results, dots show where mass 
incompleteness starts to play a role. The error bars represent bootstrap errors. The dashed and dotted lines represent results from Damen 
et al. (2009) and Zheng et al. (2007), respectively. SDSS data were used to include a local data point {triangles). Linear fits to the 
FIREWORKS results are displayed using gray dash-dotted lines. The fits show that the SSFR increases with z at a rate that appears 
independent of mass. 



the semi- analytic models of de Lucia & Blaizot (2007) 8 . 
In a qualitative way, the results of Guo & White (2008) 
concerning star formation are consistent with the ob- 
served trends in Figure O They find that the growth 
rates through star formation increase rapidly with red- 
shift, for all stellar masses. In addition, they find 
that specific growth rates through star formation are 
smaller in high-mass than in low-mass galaxies at all 
redshifts. Although qualitatively their model predictions 
agree with our results, quantitatively this is not the case. 
We determine dimensionless growth rates for our sample 
similar to Guo & White (2008). The growth rate is de- 
fined as GR sf = SFR/M* ■ t(z): the sSFR multiplied 
with the age of the universe at the redshift of observa- 
tion. We determine mean values of this growth rate at 
all redshifts, in the same mass bins Guo & White (2008) 
use. The results can be seen in Figure El The gray dots 
and black line represent the FIREWORKS growth rates 
and their binned mean. The results from Guo & White 
(2008) are shown in red. For comparison with the lo- 
cal universe we used SFRs and masses from the SDSS 
to determine growth rates at z = 0. We included both 
estimates from SED fitting (J. Brinchmann, private com- 
munication, black star) and from emission lines (Brinch- 
mann et al. 2004, magenta star). The SAMs are tuned 
to fit the local universe and, except for the highest mass 
bin, the agreement is excellent. 

8 Please note that the de Lucia &; Blaizot (2007) model assumes 
a Chabrier IMF. This will not significantly affect the comparison, 
since it greatly resembles the Kroupa IMF that we apply. 



Although the growth rates agree at redshift z = 0, at 
higher redshifts the observed growth rates increase much 
more rapidly than the model growth rates. The over- 
all offset between redshifts z = 1 and z — 2 is a fac- 
tor of ~6.3 for all masses. Per mass bin, we find off- 
sets by factors of 4.0 (8 • 10 10 M Q - 16 ■ 10 10 M Q ), 6.2 
(4 • 10 10 A/ Q - 8 • 10 10 M Q ), 6.0 (2 • 10 10 M Q - 4 • 10 10 M Q ) 
and 9.1 (10 10 M Q -2-10 10 M Q ) at redshifts z= 1-2. This 
agrees well with recent results in the literature. Elbaz et 
al. (2007) and Daddi et al. (2007a) found an excess in 
observed sSFRs with respect to those predicted by the 
Millennium Simulation of a factor of ~ 2.5 (z ~1) and ~ 
4 (z ~2), respectively, for galaxies of 3 • 10 10 M Q . Recent 
work by Santini et al. (2009) shows a similar trend where 
the Millennium Simulation undepredicts star formation 
activity by a factor 3-5. Though both these studies and 
our own are based on observations of the same field, there 
are a differences in sample selection, determination of 
SFRs and photometric redshifts. Given the wide variety 
in methods and samples used, it is remarkable that the 
observations converge to consistent results, namely high 
sSFRs with respect to model predictions. 
This discrepancy could have its origin in the derivation 
of our SFRs based on the MIPS 24 ^m fluxes. To test 
this we include the growth rates derived from the SFRs 
from our SED fits (blue lines) . The growth rates derived 
from the bolometric fluxes and the SED fits are generally 
consistent with each other, except at the highest masses. 
This result is remarkable and it argues against the hy- 
pothesis that a simple conversion error in the calculation 
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Fig. 3. — Dimensionless growth rate through star formation as a function of redshift. The different panels represent different mass bins 
in units of 10 10 Mq. The gray dots are the FIREWORKS values, the black line is their mean. Open circles with arrows denote upper 
limits set at the 2a level. Overplotted in red are the results from Guo & White (2008), based on the Millennium Simulation. Additional 
observational measurements are added, based on SFRs determined from SED fitting (blue line) and the results from SIMPLE (green line, 
Damen et al. 2009). To represent the local universe, both a photometric (black star) and a spectroscopic (magenta star) measurement 
from SDSS were included. Note the large offset between the model and the observations between redshift z=0 and z=l. The model and 
observations both show an increase of the growth rate with redshift, but for the model this trend is more gentle. 



of the SFR causes the problem. It is interesting to note 
that the SED fits generally give lower sSFRs for the high- 
est mass galaxies than the sSFRs based on UV + 24 /*m. 
AGN may be partially responsible (Daddi et al. 2007b). 
We will return to the effect of errors in the SFR deter- 
mination in the next section. 

A different cause could be that the CDFS suffers from 
cosmic variance. We added in the results from the SIM- 
PLE survey, which is based on observations on the E- 
CDFS and has ~6.5 times more area (green lines). Both 
fields show similar results. As an additional exercise, we 
redetermined the SIMPLE growth rates excluding the 
CDFS and for the CDFS only, to see whether the differ- 
ent areas show substantial differences in growth rates. 
We find that for all masses and at all redshifts, the 
growth rates for both areas are consistent with each other 
within la. Field to field variation, therefore, does not 
seem to cause a significant error. 

It is interesting to note, that although the steepness of 
the observed slope is not reproduced by the model, the 
slope is the same for galaxies of low and high mass. To 
illustrate this, we used the model growth rates to derive 
sSFRs for the Guo & White (2008) data. The results are 
shown in Figure [4] Different colors represent the same 



mass bins used in Figure [3] Both for the model and ob- 
servations, the rate of increase of sSFR does not seem to 
be a strong function of mass. 

5. DISCUSSION 

We use the FIREWORKS survey to investigate the 
star formation history in the CDFS. We find that the 
sSFR is an increasing function with redshift for all 
masses. The slope of the trend is similar for galaxies 
of all masses, although some uncertainties remain due 
to incompleteness at the low-mass end and low number 
statistics in the highest mass bin ( > 1O 11 M ). These 
findings are in agreement with results from previous 
studies (Martin et al. 2007, Zheng et al. 2007, Perez- 
Gonzalez et al. 2008, Damen et al. 2009). 
We also compare our results to model predictions of Guo 
& White (2008), who use the galaxy formation model 
of de Lucia & Blaizot (2007) based on the Millennium 
Simulation and find that for both the model and the 
observations the growth rate increases with redshift 
for all masses. However, the model fails to reproduce 
the steep rise that is observed in the data, particularly 
between z = and z = 1, a redshift range where 
observations are most reliable. The overall value of the 
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Fig. 4.— sSFR versus redshift for the Guo & White (2008) 
results. Different colors represent the mass bins used in Figure 
[3] The slope of the evolution of sSFR with redshift is nearly in- 
dependent of mass, which agrees with the FIREWORKS results. 
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Fig. 5. — Dimensionless growth rate as a function of redshift 
for galaxies with masses ranging from 2 ■ 10 10 Mq to 4 ■ 10 10 Mq. 
The black solid line shows the FIREWORKS result when a Kroupa 
(2001) IMF is applied. The dashed-dotted lines, dashed and dot- 
ted lines show the FIREWORKS growth rate based on an evolv- 
ing IMF, according to the parametrization of van Dokkum (2008; 
vD08), Dave (2008; D08), and Wilkins et al. (2008; W08), respec- 
tively. The corrected values based on the IMFs of DOS and W08 (in 
blue) are lower limits, since they do not include the effect an evolv- 
ing IMF has on the stellar mass. Introducing a time-dependent 
IMF decreases the discrepancy between the observations and the 
simulated results from Guo & White (2008; red solid line), but it 
does not completely resolve it. In particular, the steep increase in 
observed growth rate at low redshift (z = - 1) is still evident. 



growth rate is higher in the observations than in the 
models, by a factor of ~6.3 on average. This discrepancy 
between the model and the data is caused by problems 
with the observations, the model, or both. 
First of all, field to field variation may play a role. We 
have already discussed in Section [4] that this is unlikely 
to be a large effect. An additional factor of uncertainty 
resides in our estimate of the SFR based on MIPS 
fluxes. AGNs would contribute to the 24 /mi emission 
and generate overestimated SFRs. Although we removed 
all X-ray detections from our sample, dust-obscured 
AGNs that show strong 24 fim flux can not be iden- 
tified. Furthermore, systematic uncertainties in the 
conversion from MIPS fluxes to total IR luminosities 
are of the order of a factor 3 at z ~ 2 (Labbe et al. 
in preparation). In addition, we can not be sure that 
the local IR templates on which we base our conversion 
still hold in the z > 1 universe. Marcillac et al. (2006) 
suggest that local templates are still accurate out to 
z ~ 1. Papovich et al. (2007) argue that at z ~ 2, 
local templates overpredict the total IR luminosity by a 
factor of a few. However, their local templates include a 
luminosity dependent conversion. In contrast, we apply 
a single template (without luminosity dependence), 
which results in smaller systematic offsets. In addition, 
we found good agreement between the average SFRs 
based on MIPS flux and SED modeling for the low- to 
intermediate-mass galaxies. It is interesting to compare 
and investigate other measurements of the SFR. SFRs 
based on H Q , for instance, agree quite well with SFRs 
based on SED modeling (Forster Schreiber submitted). 
Furthermore, growth rates based on H Q show significant 
disagreement with model predictions (Gcnzel et al. 
2008), similar to what we find. Concluding, it is unlikely 
that the full discrepancy is caused by the systematic 
uncertainties in deriving SFRs. 

A third reason why our SFRs are large, could be the 
assumption of a non-evolving IMF. Recently, there 
have been several authors who suggest that the IMF 
is not constant, but instead evolves with redshift. 
Van Dokkum (2008; vD08) found that reconciling 
the evolution of color and mass-light-ratio of massive, 
early-type cluster galaxies favors a top-heavy IMF. 
Dave (2008; D08) argues that an evolving IMF brings 
the observed evolution in the relationship between 
SFR and stellar mass in better agreement with model 
predictions. Finally, Wilkins et al. (2008; W08) claim 
that the assumption of a simple model for an evolving 
IMF significantly reduces the discrepancy between the 
integrated SFH and stellar mass density measurements. 
We used the effect an evolving IMF would have on SFRs 
(Figure 3 from W08, Figure 3 from D08) and SSFRs 
(Figure 12a and Figure 14b from vD08) to redetermine 
our growth rates and show the results in FigureEf. It 
is important to note that for the W08 and D08 growth 
rates (shown in blue), we only included the effect an 
evolving IMF has on SFR, and not on mass. For this 
reason the blue curves should be regarded as lower 
limits. The overall factor with which the growth rate 
decreases when an evolving IMF is applied is ^1.3 for 
van Dokkum (2008), <1.8 for Wilkins et al. (2008) and 

9 The model predictions were kept to the values based on the 
non-evolving Chabrier IMF 
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< 2.4 for Dave (2008). Although this helps reducing the 
discrepancy between the instantaneous SFR and SFR 
based on the mass function (W08), the steep increase 
between z = and z = 1 is still intact. 
Finally, the use of photometric redshifts introduces 
some uncertainty to the SFR estimates. We note, 
however, that the growth rates of SIMPLE (green line 
in Figure [3]) are based on the photometric redshifts of 
COMBO-17, which are highly accurate out to z ~ 0.8. 
Since the SIMPLE results show the same strong slope 
as the FIREWORKS growth rates, we do not expect 
photometric redshift errors to be a significant cause for 
the discrepancy between the model and observations. 
None of the effects we discussed above can provide 
a clear-cut solution to explain the difference between 
the model and observations, particularly the observed 
steep rise in growth rate between z=0 and z=l. The 
discrepancy could be caused by a flaw in the models, e.g. 
concerning for instance the gas supply to the galaxy. 
The current way stars are formed in SAMs is by heating 
gas to the virial temperature after which it gradually 
cools to the temperature where star formation can 
commence. Dekel et al. (2009) show that this process is 
too slow to explain the fact that massive galaxies have 
already formed most of their stars at high redshift. They 
introduce a new way of star formation, where cold gas 
enters the dark matter halo through filaments and can 



start forming stars immediately. It would be interesting 
to investigate how this implementation of star formation 
would alter model predictions and whether it would be 
able to reproduce the steep rise of the growth rate at 
low redshifts. 

Despite the differences we identified between our obser- 
vations and the model values of Guo & White (2008), 
we also found an interesting agreement. It is remarkable 
that both observations and models find that the increase 
in sSFR is the same for low- and high-mass galaxies. 
There is no evidence for a "break" redshift at any of the 
masses studied here. Any modification in model recipes 
would have to maintain this mass independence. 
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